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Jbnsson, Yoshino and Nordblad reply: In their com- 
ment to our work 0, L. Berthier and J. P. Bouchaud 
report results of temperature-shifts AT in a numerical 
simulation of a 4 dimensional EA Ising spin-glass model 
using the same thermal protocol as in lj. Their scaling 
analysis yields apparently similar result to ours. 

In our experiments the data are limited to a narrow 
temperature region |AT| < 0.02 T g . The expected sat- 
uration of L e fi to the overlap length Tat El is not ob- 
served for these small values of |AT|. In this reply, we 
instead point out that our experimental results can be 
understood in terms of temperature chaos as partial re- 
juvenation in a weakly perturbed regime. 

Within the droplet theory 0] , it is assumed that there 
are dangerously irrelevant droplets whose free-energy 
gaps Fl are arbitrary small. The gap can be smaller than 
dF with probability p{0)dF/T{L/ L ) e , where p(0) > 
0. Consequently a perturbation which amounts to a 
small free-energy gain of order AU(L/L ) a can induce 
a droplet excitation with a probability p(L/Lau) ~ 
p(0)(T/L A c/) c . Here L AU = L (AU /T)" 1 ^ is the over- 
lap length with the chaos exponent C = ct — 9. In the 
case of T-shifts, AC/ = fc B |AT| and a = d s /2, where 
d s is the surface fractal dimension of the droplets. Im- 
portantly, the probability p(L/Lau) is non-zero in the 
weakly perturbed regime L < Lau and increases with L 
if C > smoothly connecting to the strongly perturbed 
regime L > Lau where p{L/Lau) = 10- 

Supposing that Lat > Tji^w), all length scales 
smaller than the domain size Txi(i w ) belong to the 
weakly perturbed regime. A finite region of size L within 
a domain can then become out-of-equilibrium with the 
probability piLjL at)- This reflects the breaking up of 
domains into smaller ones, yielding a broad distribution 
of domain sizes up to the maximum T max = Tj^tw)- 
Then a simple scaling argument gives that the effective 
domain size surrounding an arbitrary point of the sys- 
tem on average reaches, L e g = L at F(L m&x )/ Lat) with 
F(x) — x — cx 1+< ^. Here, the 2nd term in F(x), with c be- 
ing a positive constant yields a correction term, which is 
the first order of p(0) AC/ /T to the no-chaos limit of L e ff ■ 
Higher order terms can be neglected for small enough 
AT|. It can be seen in Fig. ^that the ansatz fits our 
data for \AT\ < 0.4 K. 

In the simulations of Berthier-Bouchaud 0, length 
scales L < 5 are probed while Lat ^20. In the 
spatial correlation function is found to decay to only af- 
ter LTiit-w) being consistent with L max = Tj^w)- Then 
the rejuvenation found in can be understood as par- 
tial rejuvenation due to the temperature chaos effect in a 
weakly perturbed regime. We note that emergence of the 
temperature chaos effect has been confirmed in the same 
model as in with modest system sizes @- Indeed we 


have found that the ansatz proposed above also fits the 
result of Berthier-Bouchaud very well with c = 0.15 
and ( = 1. However, the largest T-shifts used in in- 
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FIG. 1: Scaling plot of L efi with L AT = L {a\AT\/ J)~ 1/z 
with a = 5 and 1/C = 2.6. The fit is due to the scaling ansatz 
for the weakly perturbed regime explained in the text. The 
open/filled symbols are for positive/negative T-shifts. 


volve temperatures in the range 0.8 — 0.9T g , within which 
numerical (but not experimental) length/time scales are 
dominated by critical fluctuations as shown by Fig 1. of 

We conclude that the results reported in [l| can be con- 
sistently understood in terms of the temperature chaos 
effect in a weakly perturbed regime. We are not aware of 
any other concrete scenarios which predict the symmet- 
rical effect of positive and negative T-shifts reported in 
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